A Matlab-based finite-difference solver for the Poisson problem with mixed Dirichlet-Neumann boundary conditions
نویسندگان
چکیده
A Matlab-based finite-difference numerical solver for the Poisson equation for a rectangle and a disk in two dimensions, and a spherical domain in three dimensions, is presented. The solver is optimized for handling an arbitrary combination of Dirichlet and Neumann boundary conditions, and allows for full user control of mesh refinement. The solver routines utilize effective and parallelized sparse vector and matrix operations. Computations exhibit high speeds, numerical stability with respect to mesh size and mesh refinement, and acceptable error values even on desktop computers. PROGRAM SUMMARY/NEW VERSION PROGRAM SUMMARY Manuscript Title: A Matlab-based Finite-difference Solver for the Poisson Problem in Two and Three Dimensions, Optimized for Strongly Heterogeneous Boundary Conditions Authors: Ashton S. Reimer and Alexei F. Cheviakov Program Title: FDMRP 1.0 Finite-Difference with Mesh Refinement Poisson equation solver Journal Reference: Catalogue identifier: Licensing provisions: GNU General Public License v3.0 Programming language: Matlab 2010a Computer: PC, Macintosh Operating system: Windows, OSX, Linux RAM: 8GB (8,589,934,592 bytes) Number of processors used: Single Processor, Dual Core
منابع مشابه
A fast spectral/difference method without pole conditions for Poisson-type equations in cylindrical and spherical geometries
A simple and efficient FFT-based fast direct solver for Poisson-type equations on 3D cylindrical and spherical geometries is presented. The solver relies on the truncated Fourier series expansion, where the differential equations of Fourier coefficients are solved using second-order finite difference discretizations without pole conditions. Three different boundary conditions (Dirichlet, Neuman...
متن کاملMaximum-norm stability of the finite element Ritz projection under mixed boundary conditions
As a model of the second order elliptic equation with non-trivial boundary conditions, we consider the Laplace equation with mixed Dirichlet and Neumann boundary conditions on convex polygonal domains. Our goal is to establish that finite element discrete harmonic functions with mixed Dirichlet and Neumann boundary conditions satisfy a weak (Agmon-Miranda) discrete maximum principle, and then p...
متن کاملNumerical resolution of large deflections in cantilever beams by Bernstein spectral method and a convolution quadrature.
The mathematical modeling of the large deflections for the cantilever beams leads to a nonlinear differential equation with the mixed boundary conditions. Different numerical methods have been implemented by various authors for such problems. In this paper, two novel numerical techniques are investigated for the numerical simulation of the problem. The first is based on a spectral method utiliz...
متن کاملA Consistent and Accurate Numerical Method for Approximate Numerical Solution of Two Point Boundary Value Problems
In this article we have proposed an accurate finite difference method for approximate numerical solution of second order boundary value problem with Dirichlet boundary conditions. There are numerous numerical methods for solving these boundary value problems. Some these methods are more efficient and accurate than others with some advantages and disadvantages. The results in experiment on model...
متن کاملImproving the Rate of Convergence of ‘high Order Finite Elements’ on Polyhedra I: a Priori Estimates
Let Tk be a sequence of triangulations of a polyhedron Ω ⊂ Rn and let Sk be the associated finite element space of continuous, piecewise polynomials of degree m. Let uk ∈ Sk be the finite element approximation of the solution u of a second order, strongly elliptic system Pu = f with zero Dirichlet boundary conditions. We show that a weak approximation property of the sequence Sk ensures optimal...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- Computer Physics Communications
دوره 184 شماره
صفحات -
تاریخ انتشار 2013